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ON HERBICIDAL ACTIONS 


I. OPATOWSKI AND ALICE M. CHRISTIANSEN 
UNIVERSITY OF MICHIGAN AND WRIGHT JUNIOR COLLEGE, CHICAGO, ILLINOIS 


A theoretical interpretation of the herbicidal action of the 2,4-Dichlo- 
rophenoxyacetic acid is suggested. 


It was the idea of P. Jordan (1941) to interpret certain biologi- 
cal effects as due to a random event consisting of a penetration of a 
molecule of an external agent into a vital center of the organism (v. 
Schelling, 1942). It is the purpose of the present note to point out 
that the experimental results of J. van Overbeek, L. E. Gregory and 
I. Vélez (1946) on the herbicidal action of water solution of 2,4-Di- 
chlorophenoxyacetic acid are in agreement with these concepts, if 
such action is attributed to a single molecule or to a single aggregate 
of molecules of the acid reaching a particular center of the plant. No 
direct experimental evidence of such mechanism can be obtained, at 
least with the presently known techniques. Some indication of its 
possibility may be seen in the known general fact of translocation of 
the chlorophenoxy compounds within the plants (Zimmerman and 
Hitchcock, 1942). Other indications are found in the specific experi- 
ment of J. van Overbeek, L. E. Gregory, and I. Vélez (1946) consist- 
ing of painting a leaf of shoots of commelina with a solution of 2,4-D, 
which produced the death of a stem region 10 cm below the treated 
leaf and did not affect the stem above the leaf. According to these 
authors, 2,4-D does not burn the foliage but acts internally on a well- 
defined region of the plant. 

The following theory is suggested for an interpretation of the 
quantitative results of J. van Overbeek and his co-workers: Let a 
homogeneous aggregate of plants be sprayed with a solution of con- 
centration C, the concentration being measured in volume. C is the 
total volume of the molecules of the acid per 1 cm®* of the solution. 
Let S be the total volume of the above mentioned sensitive vital cen- 
ters of the plant and R the volume of a region in the plant contain- 
ing those centers and such that the solution does not penetrate in the 
plant beyond that region. The region R may be the whole volume of 
the plant. Let V be the volume of the solution in R. If M cm? is the 
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volume of one molecule of the acid, there are C/M molecules of the 
acid per 1 cm of the solution, i.e., the average volume of the solution 
containing 1 molecule of the acid is M/C cm’. If V = R, i.e., if the 
solution penetrates the whole region R and consequently the whole 
volume S, then SC/M is the average number of molecules of the acid 
in S em’. If the volume V of the solution in R is smaller than R, then 
in S cm? of R there is an average of (V/R)S cm’ of solution and the 
average number of molecules of the acid in S cm’ is VSC/(MR). The 
maximum possible number of molecules of the acid in S cm’ is S/M. 
Consequently, if we divide the volume S into S/M parts of M cm* 
each, the probability that a molecule of the acid reach one such par- 
ticular part is represented by 


VSCM> R*/SM1=VC/R. (1) 


The probability that no molecule of the acid reach that particular part 
of S is 1—VCR-, and the probability that no acid molecule reach any 
of the above parts of S is 


(1— VCR-)8/™ , (2) 
The probability P that at least one molecule of the acid reach S is: 
P=1— (1—VCR")§”™., (3) 


Since (1+y)* = exp(yz), for small y and large z, we obtain under 
the assumptions of small concentration C and a large sensitive volume 
with respect to the volume of one acid molecule: 


P =~1—exp(—VCSR" M-), (4) 


The above reasoning follows to some extent the pattern used by 
M. A. Lauffer and W. C. Price (1945) in the study of infection by 
viruses. The theory presented here disregards the possibility that 
the solution may be eliminated from the plant through natural pro- 
cesses and that new solution may take its place. The possibility of 
recovery of the sensitive volume to its original condition is also not 
taken into account. The probability P represents what J. van Over- 
beek and I. Vélez (1946a) call the “mortality” of the plant, if the 
time between the occurrence of the critical event (e.g., combination 
of a molecule of the acid with a particular organic molecule of a 
plant cell) and the appearance of the external effect (e.g., destruction 
of a growing region of the plant) is smaller than the interval of time 
between the spraying of the plants and the counting of the dead ones, 
and if the latter interval of time is sufficiently long so that no fur- 
ther translocation of the solution within the plant occurs. 
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A comparison of the values of P calculated from the above for- 
mula with the experimental results of J. van Overbeek and I. Vélez 
(1946a) is given in the table below. Those authors published their 
results in form of smooth thick curves which give for each value of 
concentration a range of values of mortality. The curves are clearly 
of an exponential type, as it is required by the probabilistic theory 
based on a “single effective event” type of process. It is seen from 
the table that the fit is good, although working directly on experimen- 
tal figures rather than on the smoothed curves could lead to more 
definite conclusions. 


In the following table C denotes the concentration of the acid and 
he) er SAIS 


COMMELINA AMARANTHUS CYPERUS 
S oS S R 3 s 
=? zo) 
Ba, 28 ga, $8 Ba, 25 
2 Ee eal Ee a =. ol 
aS) Lx eo) 02 a) SES 
@ & @ &9 — &0 Ss 
S Fe és Be és pe 6 | OS 
S Be ae race Le ape aE 
0.25 | 0.80-0.87| 080 | 047-050] 0.48 | 018-027] 0.28 
0.375 | 0.90-0.92] 0.91 | 0.60-0.63 | 0.63 | 0.29-0.32| 0.32 
0.50 0.94-0.96] 0.96 | 0.68-0.73 | 0.73 | 0.39-0.41 | 0.40 
0.75 0.98-0.99|; 0.99 | 0.84-0.87| 0.86 | 0.53-0.56 | 0.54 
1.00 0.91-0.93 | 0.93 | 0.64—0.67 |- 0.65 
1.25 0.94-0.97 | 0.96 | 0.72-0.74 | 0.78 
1.50 0.98-0.99 | 0.98 | 0.79-0.81 | 0.79 
2.00 0.88-0.90 | 0.88 
2.50 0.93-0.96 | 0.98 
3.00 0.96-0.98 | 0.96 


The knowledge of k makes it possible to estimate the total sen- 
sitive volume S in a plant: 


S=(R/V)Mk or S2Mk. (5) 
If the experiments are carried out in condition of saturation, that is, 
with such an amount of solution that any further increase of this 
amount does not modify the mortality, then V = R and S= Mk. The 
sensitive volume is then simply k times the volume of one molecule of 
the acid. If this condition was fulfilled in the experiments in ques- 
tion, then the sensitive volume in each plant was from about 1000 to 
about 6000 times the size of a single molecule of the 2,4-Dichloro- 
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phenoxyacetic acid: However, an estimate of the volume of a mole- 
cule is a difficult problem. Moreover, a definition of that volume in- 
dependently of the nature of the sensitive center does not seem pos- 
sible, because both the volume of the molecule of the acid and that of 
a molecule of the sensitive substance, as they are used in the present 
theory, are regions of space bound to those molecules and charac- 
terized by the fact that a mutual contact of these two regions leads 
to the event observed, for instance, through a single chemical reac- 
tion or through a chain of such reactions. Such definition may be 
still incomplete, because it does not take into account the fact that 
the occurrence of a certain chemical reaction may depend not only on 
an overall distance, but also on a particular reciprocal orientation of 
the two molecules. 

In any way, routine computations of the molecular volume (cf. 
Glasstone, 1940) give a figure of about 300 A* per molecule of the 
2,4-Dichlorophenoxyacetic acid. Very rough estimates based on in- 
teratomic distances (cf. Stuart, 1934) give a round figure of 50 A®*. 
If 300 A* is taken for M , the total sensitive volume existing in each 
plant is S 2 300k A’, which is about 8 X 10-19 em? to 2 * 10°78 cm, 
depending on the species of the plant. This is an extremely small 
amount if compared, for instance, with the 10°, gr. of auxin existing 
in yeast per each gram of dry weight (Kogl and Haagen Smit, 1931; 
Kogl, 1933). 

Dr. H. D. Landahl has kindly pointed out to the writers that the 
above theory does not take into account the existence of a threshold 
concentration of 2,4-D. Such concentration has been actually given 
by P. W. Zimmerman (1943). In this connection it should be kept in 
mind that the threshold concentrations are extremely low. Those giv- 
en by P. W. Zimmerman are some six thousand times smaller than 
the maximum concentration used by the same author in the same set 
of experiments. Experimentation with such low concentrations is 
very inaccurate unless a very large sample is used. 

The authors are indebted to Professor H. D. Landahl of the Uni- 
versity of Chicago for his enlightening criticism, to Professor F. G. 
Gustafson of the Department of Botany of the University of Michi- 
gan for suggesting invaluable references and to Professor K. Fajans 
for his kindness in discussing some general questions of chemistry 
connected with the present subject. 
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ON MATRICES OF NEURAL NETS 


TELSON WEI 
ACADEMIA SINICA 


The structure of a complete or incomplete neural net is represented 
here by several matrices. The activity equation of the net follows in a 
general form. A chain or cycle is defined as a neural structure whose 
connection matrix is unitary. We can compute the number of simple 
chains by a recurrent formula. 


In a paper by H. D. Landahl and Richard Runge (1946), the ac- 
tivity of a neural net is represented in terms of a matrix vector equa- 
tion in which the matrix represents the structure of the net. Here we 
represent the structure of a net by several matrices which have their 
own peculiar significances. 

We shall make the following physical assumptions for our neuron- 
nets: 


1) The activity of the neuron is of the all-or-none character. 


2). A neuron can be excited into all-activity either by the actions 
of other neurons connected to it or by an external stimulus. An 
afferent neuron is one stimulated by an appropriate external stimu- 
lus. An internal neuron is stimulated only by other neurons. An 
internal neuron fires if and only if a sufficient number of end bulbs 
terminating on it fire within the period of latent addition. The 
minimum number of terminal bulbs which must fire within the 
period of latent addition in order to make the neuron fire is called 
the threshold of the neuron. The threshold is considered as a con- 
stant independent of the previous activity of the neuron. The all- 
activity of a neuron is always of the same intensity through every 
branch of its axone. 


3) The all-activity of a neuron lasts throughout its time of action. 


A) The synaptic delay is the same for all neurons in a neuron-net 
and is the only significant delay in the nervous system when it 
is compared with the period of latent addition, the time of action 
and the conduction time. For convenience, the interval of synaptic 
delay will be considered as a unit time. 


5) No neuron acts upon itself. 
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6) In the configuration of our network we think of the resultant 
effect of all excitatory and inhibitory synapses from neuron Nj; 
directly to neuron N; as a single action through a single fiber, ||t j\ |. 
The action is positive, negative or zero according to whether it is 
purely excitatory, purely inhibitory or has no effect at all. Any 
action from the same source to the same object directly, i.e., no neu- 
rons lying in between, is also considered as through a single fiber. 


General Activity Equation. By a complete net we mean that no 
other neurons or other fibers outside of this net are associated with 
it; otherwise, it is incomplete. Let N be a net of m ordered neurons 
(N,, N2,-::, Nn), complete or incomplete. Let 


s,(t) = the intensity of external stimulus acting upon N; at the 
time t, 


ni; = the non-negative intensity of stimulus through ||7 7|| , 


+1 ais 
dij=}—1 when |iij||;<0, 
0 As, 
ti; = di; ni; , 
h; = threshold of N; , 
According to our assumption, di; = ni; = ti; = 0 fori=1, 2 preg hs 
We define 


a(t) = (a, (t), M(t), ---, @n(t)) 


=the activity vector of N, as defined by Landahl and 
Runge, 


E = (nj;) =the intensity matrix of N , 

D = (d;;) = the connection matrix of N , 
T= (t;;) = the structural matrix of N , 
1h eee 0 | 


= the diagonal inverse thresh- 
old-matrix of N. 


The scalar quantity 


1 jen : 
v;(t) as (si (€) + ¥ a;(¢—1) nj; di) 
a gut 
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measures the excitation state of N; at time t. The vector 
v(t) = (v,(t),-+-, vn(t)) = (a(t—1) T+ 8(t))H 


measures the excitation pattern of N, where s(t) = (si (4), +++, Sn (¢)). 
A quantity q can sometimes be normalized in the following way: 


{qe when "921; 
[a =| j when qg<1l. 


A matrix is said to be normalized and denoted by the similar no- 
tation if we normalize all its elements. 
We have thus the activity equation of N, 


a(t) = [v(t)] =[(a(t—1)T + s(t)) A]. (1) 


Chains and Cycles. A complex C,, of m neurons and m associated 
fibers in a net constitutes a simple (structural) chain of order m if 
each neuron receives fibers from only one neuron and sends fibers to 
only one neuron. A simple chain is called a cycle if it cannot be re- 
duced to two or more simple chains. 

A mathematical representation of simple chains and cycles can 
be obtained from the following theorems: , 

THEOREM 1. The necessary and sufficient condition for a com- 
plex C,, in a net to be a simple chain is that the connection matrix D 
of Cm is unitary, i.e., DD* = D*D = e, where D* is the transpose of 
D and e is the unit matrix of orderm. (See Pontrjagin 1939, p. 108.) 

In fact, Cn is a simple chain if and only if there is one and only 
one element different from zero, viz., +1 or —1, in each row and in 
each column of its connection matrix D. Therefore, if and only if 


Dd di; dix = bjx (6;=—1, 6i;=0 fori1#7). 
4=1 
THEOREM 2. A simple chain of order m can be decomposed 
into a finite system of disconnected cycles of which the sum of orders 
is equal to m if and only if its connection matrix is reducible. Every 
cycle of order r is represented by an irreducible unitary connection 
matrix of orderr. 
By rearranging the order of neurons in a simple chain of order 
m , we can always have the connection matrix in the new form 


~ oy Fads O 
t Dt=(4 “Jk 


where t is non-singular. The submatrices a and b are also unitary 
and represent two disconnected simple chains. The sum of the orders 
of aand bis m. If aor b is still reducible, it in turn may be further 
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decomposed. In this way we get finally a finite system of irreducible 
submatrices, each of which is unitary and represents a cycle. 

In order to evaluate the number of simple chains and cycles in 
a net, we define the diagonal submatrices of a matrix as those formed 
by taking some rows and the corresponding columns from it. From 
Theorem 1 follows immediately 


THEOREM 3. The number of simple chains including all cycles 
in a net is equal to the number of unitary diagonal submatrices of the 
connection matrix of the net. 


In practical calculation for some special cases we formulate the 
simple chains and cycles in another way. Regarding ||i j|| as 7 sub- 
stituted by 7, we get a two-row notation of a substitution or permu- 
4s ley ba 
Jao Jao°e*s Im 
stituted by 7.. We often employ the one-row notation of cycles so 
that each substitution can be decomposed into several disconnected 
cycles. (See Dickson 1930, p. 139). Conversely every substitution 
of order m of subscript-pairs ||i 7|| # 0 gives a simple chain of order 
m. 


tation for a simple chain of order m ( ) , where 7, is sub- 


LEMMA 1. Each non-zero term of the determinant of a diagonal 
m-row square submatrix of D gives a simple chain of order minN. 

Because in every non-zero term of such a determinent D,, , which 
is the product of m factors, +1 or —1 formed by taking one and only 
one element from each row and each column of D,, , all subscript-pairs 
of these factors form a substitution. 


LEMMA 2. The determinant |D| has at most 
a, = (n—1) ((n—1) !— e,4"°) (2) 


non-zero terms, where ex* is determined by successive deduction from 
the formula 


Cn* = (m— 1)!+ (k— 1) Cm’? + (m— k) Cm (3) 


fork Smandk > 2 with en® =0, ent = (m—1)!, en? = (m—1)! + 
(m — 2)e,,. 

We consider first a determinant D,, of order m in which k diago- 
nal elements are zero and all other elements are non-zero. Let e,* be 
the number of zero-terms in the expansion of Dn. By successive in- 
terchanges of rows and columns we see first that e,* is independent 
of the distribution of those zeros in the diagonal line and that ‘any 
determinant having but k& zero-elements, of which no two lie in the 
same row or the same column, has exactly e,,* zero-terms in its ex- 
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pansion. Formula (3) is thus easily derived by induction from Dy» 
with zero-elements in simple arrangement. Consequently the determi- 
nant |D| has non-zero terms in a number not more than 


ad, =n! — e,"=n!— (n—1)!— (n—1)@,,."° 
ma ag LAC i Lat Ong a) 


From these lemmas follows immediately 


THEOREM 4. The net N has at most a, simple chains of order 
n, where a, is determined by (2). 
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A MATHEMATICAL STUDY OF ORGANISM GROWTH AS AN 
APPROACH TO THE CANCER PROBLEM 


RUSSELL H. KESSELMAN 
THE JEFFERSON MEDICAL COLLEGE OF PHILADELPHIA 


_ A mathematical expression is described giving the growth rate of a 
tissue region as a function of the determining biological factors. It is 
observed that within the limitations of prediction of this expression, a 
change of growth to the neoplastic type can be expected under certain 
conditions of the variables. The growth equations used here are dif- 
ferent from others usually mentioned in that they attend to both the 
positive and the negative growth phases occurring in the course of or- 
ganism development. The expression derived enables the prediction of 
an increase in incidence of cancer with age, predicts a decrease in the 
rate of tumor growth with increasing age, and permits an hypothesis to 
be made as to the nature of the action of carcinogenic agencies. 


The study of growth may be approached on a number of different 
levels as far as environmental (both internal and external) influences 
upon the hereditary substrate are concerned. The most primitive of 
these environments concerns itself with growth and differentiation 
(Little, 1947). This is the environment which our mathematical treat- 
ment will consider. Another environmental phase of great impor- 
tance in the cancer problem, which might also be mentioned, is that 
of the phenomena controlled primarily by horomones. This study of 
the environment concerned with growth and differentiation will be 
limited to the period from the time of birth to the death of the or- 
ganism as this is seen in the higher animals. 

In this exposition cell number will be used as the parameter rep- 
resenting growth. There is considerable justification for this (Shock 
and Morales, 1942). The idea of growth is certainly very closely 
bound to the tumor problem. In fact, J. Ewing (1940) states, “A tu- 
mor is an autonomous new growth of tissue.” Speaking more directly, 
the use of cell number here as a growth parameter appears natural, 
since microscopically rapidity of growth in tumors is shown by evi- 
dence of cell division (Barnard and Robb-Smith, 1945). It should be 
understood, of course, that an increased rate of cell division is not 
actually the primary factor in tumor formation, for if it were (Berill, 
1943), as cell division continued tumor cells should get progressively 
smaller and disappear. It appears that the increased rate of cell divi- 
sion is a device to insure mechanical stability as the bulk of proto- 
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plasm is increased. Using certain assumptions, N. Rashevsky (1938) 
showed that from physical considerations a cell above a certain size 
is mechanically unstable. It might be pertinent to mention here that 
the complete solution of the growth problem demands consideration 
of changes in: (1) cell size, (2) cell number, and (3) intercellular 
substance. The most readily observed and measured change in tumors 
is that of cell number. Our primary purpose here is not a complete 
solution of the growth problem. Our interest is centered in that phase 
of growth which will best enable us to study tumors. 

In this discussion a tissue region will be defined as the whole of 
those cells of the same histological type in continuity with each other 
in the same organ (or organs). The total number of tissue regions in 
the entire body will be designated as m tissue regions and the typical 
tissue region will be known as the ith tissue region. The number of 
cells in the ith tissue will be denoted by 7; . 

Let n = the total number of cells in the entire organism. 

Then we obtain 


HEaSye. (1) 


4=1 


We shall define a constant g which will be used as a measure of 
the hereditary constitution of an organism insofar as this constitu- 
tion influences change in cell number occurring in the growth process. 
Then the amount of g in each cell of a given individual organism is 
the same as that in every other cell in that organism, since the spe- 
cific hereditary materials of the individual are possessed equally by 
all of the cells of that individual (Newman, 1932). 

It is known that when groups of cells of various tissue types are 
placed in tissue culture medium with all exogeneous influences re- 
moved, the most significant distinguishing characteristic between 
these types is the difference in their rates of cell division (Fischer, 
1946). This is true even though the characteristics of these cells from 
the point of view of genetics may have been the same, i.e., each group 
of cells came from the same organism. Therefore, we introduce a 
quantity a; which will bear a constant value in any particular tissue 
region 7. The numerical value of a; will bear an inverse relation to 
the degree of differentiation which the ith tissue region has under- 


gone. Then a; will describe the influence of differentiation on cell 
division rates. 


Let 
Gi=ag, (2) 
for a single cell of type 7. 
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The second basic quantity to be discussed is related to the con- 
cept of the individuation field in experimental embryology. It must 
be remembered that most of our knowledge about the individuation 
field deals with what in the higher animals is the prenatal stage of 
growth, but that in this particular discussion, however, we limit our- 
Selves to the post-natal stage. The existence of this field is a conse 
quence of cellular aggregation in metazoans. Body development is not 
a product of the cell but of other much more general formative forces. 
In tissue culture it has been shown that changes which the cell under- 
goes from maturity to death—growth, differentiation, and function— 
are not changes peculiar to it, but merely a reaction which it may 
suffer under the influence of different conditions in its environment 
(Burrows, Jorstad, and Ernst, 1926). A number of examples of in- 
teractive effects are: (1) competition of cells for nutrient substances, 
(2) mechanical effects (compression), (3) spatial effects, such as dis- 
tortion of nutritive diffusion fields of a cell by its neighbors, (4) pro- 
duction of harmful metabolites, (5) growth-stimulating agents, (6) 
the influence of specialized functions, (7) organization, and (8) ac- 
cumulation of products which are not dissolved out under the condi- 
tion of organization found normally in the body of the organism. We 
will use the symbol J; to represent the effect of such interactive influ- 
ences on cell division rates in the 7th tissue region. 

The quantities G; and J; will be assumed to have additive prop- 
erties so that we can let 


dt mG,—Ji, (3) 
where ¢ is taken as time and where G; is multiplied by 7:, since we 
are dealing with a tissue region of n; cells. 

The quantity J; is a function of time, and by the accumulation of 
J;, when cells aggregate in sufficiently large numbers, their growth 
is gradually slowed down and may eventually be given a negative 
value. The quantity J; is usually an increasing function of age. Ana- 
tomically such increase might be said to be observable in the decrease 
in cellular elements and with the increase in intercellular substance 
coming with age. Changes in tissue fluids with aging of the metazoan 
body are inevitable. They are probably of many kinds and result 
sooner or later in alterations in the cells subjected to them. The cells 
fail to function as they should and the whole cellular community may 
be invalidated or die as a result (Cowdry, 1942). The same idea is 
made manifest by the work of A. Carrel and M. T. Burrows (1911), 
who found that plasma was not an optimum medium for the growth 
of cells planted in it. When diluted by two volumes of distilled water 
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to three volumes of plasma, much better results were secured. In our 
language this would be a decrease of J;. In tissue culture where en- 
vironment is frequently changed, there is no opportunity for J; to 
increase (accumulate), and cells can and do grow indefinitely. Cells 
in tissue culture demonstrate very many of the properties of neo- 
plasia. In the case of unicellular organisms and in metazoans of small 
size like hydra and planaria, gross aggregation of cells with resultant 
irreversible increase (accumulation) of J; is not present, and such 
organisms achieve a kind of immortality (Heilbrunn, 1943). 

A certain minimum value of J; (or interactive effect among cells) 
seems necessary for any growth to occur at all since, when individual 
cells or small groups of cells are brought into media where nutritive 
requirements are satisfied, growth will not occur unless a certain mini- 
mum number of cells is present (Fischer, 1946). 

Let »; be the minimum number of cells of type i which must be 
together in aggregation for any cell division to occur at all. When 
n; < »;, no cell division occurs. When n; = »;, dni/dt = 0. When 
n; >v,,dn;/dt >0. 

We will assume that 

Ji=f(n,n;,t). (4) 


We are making J; an explicit function of t, which is a departure 
from previous equations of this type (Rashevsky, 1945). In so doing, 
we recall our previous remarks on the changes of intercellular sub- 
stance and of tissue fluid occurring with age. In the final integrated 
equation which gives n; as a function of t, this will permit n; to 
reach a maximum and then decrease as the process of senescence sets 
in. In other equations which have been used heretofore, n; approaches 
a limiting maximum value with age. N. Rashevsky has pointed out 
the objectionability of the anthropopsychic concept inherent in the 
use of time in this way, which assumes that the organism knows its 
own time or carries its own clock. More desirable than the use of time 
in our equations might be the use of change in entropy due to physi- 
ological activity of the organism under discussion. Assuming J; to be 
an increasing function of entropy would result in some interesting 
theoretical possibilities. The increase in entropy accompanying the 
course of physical experience is only probable and never certain. In 
the statistical interpretation of the second law of thermodynamics 
there is always a probability that any given state will occur if one 
waits long enough. We are not prepared to say at present whether 
this is wholly comparable to the device to be used in ‘equation (43). 

Since each cell of the organism interacts with each cell of tissue 
region 7 by an average amount, we can put J; = ntf’ (1). 
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eG; Reasoning that the J; for n; cells is approximately n; times what 
it is for one cell when 7; is large, and remembering that when 7; = »; , 
dn;,/dt = 0 and when n; > »; , dn;/dt > 0 , we let 


J, = bitn(n; — »;) + Gi, (5) 


introducing the constant b; defined for the ith tissue region. Tenta- 
tively, the simplest biological concept which may be attached to b; is 
that of organization. Groups of cells of similar structure and per- 
forming similar specialized functions are called tissues. A further 
specialization or organization of the various tissues into organs is 
met within the higher forms. Each organ consists of several kinds of 
tissues built into a structural and functional unit (Carlson and John- 
son, 1941). 
Substituting equation (5) into (3) we find: 


dn,/dt = niGi wr b;tn(n; == vi) = viGi : (6) 


In a tissue culture n = 7; and b; would be far smaller in numeri- 
cal value than in the organized body. In tissue culture we let 6; = bit. 
Therefore, 

dn;/dt = 1,G; — bint (n; — 1) — iG; . (7) 


By adding up the growth rates of individual tissue regions we may 
obtain the growth rate of the body as a whole. Define the quantities 
G,b,»,andn by 


pee > bi 
G=—, b=—, (8) 
m m 
=> i, and n= Da (9) 


Then from equations (1) and (6) we obtain: 
dn/dt =nG — btn(n—v) —7G. (10) 


Unfortunately, there appears to be no other way except by using a 
power series to obtain an exact solution of equation (10). However, 
a solution of (10) is readily available if we simplify it to the form 


dn/dt = nG — btn?. (11) 


This may be at least partially justified by noting that n >> » while 
only n; > %. 

The choice of a point for making t = 0 poses a problem. Here we 
choose birth as this point. Our reasons are perhaps not totally ade- 
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quate, but include: (1) the peculiar intimacy of contact of foetal (and 
embryonic) tissues with the exterior through the medium of the pla- 
centa, whose permeability characteristics as a stimulus to growth 
might be such as to drain off accumulated J; , the destruction of this 
relationship producing a biological discontinuity at birth; (2) the 
small size of the organism during its prenatal period with resultant 
decreased irreversible increase of J;; (3) the use of conception as a 
zero point for time would make little practical difference, since early 
in the organism’s existence J; is small; (4) it is only with difficulty 
that the present theory can be used for prenatal development, since 
in such development the a;’s for the different tissue regions would 
certainly be variables rather than constants, as we have assumed 
them to be in post-natal life; and (5) prenatal life-span is very small 
compared to post-natal life-span. 

In general it should be observed that equation (11) is less de- 
scriptive of reality in the first few years of life than at other times. 
Assuming that at birth t = 0, equation (11) would give a very high 
dn/dt (or growth rate) at the birth of the individual, with dn/dt at 
birth as great as that seen in tissue culture growth. It is known, how- 
ever, that growth is very rapid in the early years of life. 

The solution of equation (11) is found to be 


n= G1 (GE 1) sce" (12) 


where ¢ is an integration constant. 


Taking the partial derivative of n with respect to ¢ in equation 
(12), we find: 


on/ot = — [G?b (Gt — 1) + ce-**]2[bG — cG-*] , (13) 
bG — cGe-¢'m = 0 for the maximum 7 , (14) 
log (¢/b) 
ae Sa RTT 
G (15) 
The maximum number of body cells is then 
Mm = G?/(b log [c/b]). (16) 


The number of cells in the body at birth is found by making t = 0 in 
equation (12). This gives: 
My = G?/(ec— bd). (17) 


The ¢ given in equation (15) is probably about 20 years, since 
this is the age when maturity is said to be first reached (Scammon, 
1942). 
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An interesting surmise suggested by equation (12) is that at the 
age which might be the absolute limit of the possible human life-span 
under optimal circumstances, there are about the same number of cells 
in the body as at birth. In senescence there are known to be many 
similarities to a youthful condition. L. V. Heilbrunn (1943) gives 110 
years as about the longest authenticated life-span. We will take 150 
years as being the limit of the human life-span (Bogomolets, 1946). 

When 


t=20, n=3.55X 10" (Arey, 1946). (18) 
When 
b= 20 , sdn/dt=—0; (19) 
Nesp — Nizises (20) 
dn/dt = Gn — btn’; (21) 
0= GB:55 X 10") —20b(3.55 X 10*)?; (22) 
G/b=7.10 X 10%. (23) 
From equations (15) and (19), we find: 
20G = log (c/b), (24) 
and from equations (12) and (20), we find: 
—) C= 150DG— ) = cere”. (25) 
From equation (24) we get: 
b/c = e%F , (26) 
Equation (25) gives us 
1= 150 ve Sn (27) 


Using equation (26) we obtain: 

VS 5 0G ee ere (28) 
Multiplying both sides of this equation by e?°¢, we find: 

e206 — e-180¢ — 150G. (29) 


Solving equation (29) graphically, we find: 
G= 0.158. (30) 


From equation (23), we obtain: 
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b= G/(7.10 X 10"), (31) 
b = 0.0223 X 10" . (32) 

Let 
=yx10™. (33) 


Introducing equation (18) in (12), we find: 


3.55 X 10% = (0.158) 2/[0.0223 X 10-4 (0.158 X 20 — 1) 
+ (y x Ee) (ene 


34 
y = 0.522 ey 


¢=—0.522-< 10-™*. 
Because of simplifying assumptions, the numerical values given 
G,b, and ¢ are orders of magnitude only. 


Using the values obtained for G, b, and c in equation (12), we 
obtain: 


(0.158) ? 


Th O88 SCA Ill (C1 Eon a oo eee (35) 
0.0223 X 10-14 (0.158¢ — 1) + 0.522 K 10-14e-0-158¢ 


The graph of equation (35) appears in Figure 1. 


Nn Eells) 


n= GY/[bGt-N+ce-**} 


3x10 


2x10" 


1x lo! 


60 
FIGURE 1 


80 in t (years) 


Remembering equation (2), we will define 


a=G/g. 
Then ig . (36) 
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G=ag. (37) 


Substituting equation (2) into (7) and equation (37) into (12), 
and then substituting equation (12) into (7), we obtain 


dn; ins = b;ta?g? (n; —= vi) 
Peas 6 blag ae cet ee Sou 
dn; b;ta*g? 
<. 9 Ste meee icsesapprieean tae ea: 
dt b(agt —1) + ce- 


b,ta’g?y; 89) 


ide iis ewe oe 
Equation (39) is a linear differential equation of the first order with 
the following solution: 


Nn; = ex eal weg! d 
; im b(agt —1) + eae ae | é 


fl an fl b,a2g°t dt 
ek dF 
b(agt —1) + ce Yi 


b;ta?g?y; 
b(agt—1) + ce 


b;ta?g? 
+evexp{— f[ “7 ___ ag |at } 
b(agt —1) + ce 


with c; the integration constant. Since the integrals in equation (40) 
can not be evaluated readily, a graphical integration will probably be 
necessary. 

Qualitatively from the conditions of the derivation of equation 
(40) we can take a guess at the form of the graph of n; plotted against 
t. There will be an increase in youth and then a decline of 7; will 
set in with increasing age. However, in certain regions, dependent on 
the values of a; , b; , and »;, a release of an adequate number of cells 
from intercellular controls will be effected at some age with result- 
ant sudden cellular proliferation occurring which will fit the descrip- 
tion of neoplastic growth. 

The cancer problem is being attacked today by a multiplicity of 
methods which are rather diverse in approach and viewpoint. We 
will make an attempt here to show how some of these may be fitted 


into our theoretical study. 


(40) 
— Vii | dt 
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It is known that in the process of growth and cell division a very 
small number of cells will occur, say two or three in tens of thous- 
ands,{which will have a heredity different from their ancestors. We 
are not speaking here of the small, continuous, accidental variations 
that are bound to occur even in the most homogeneous population, 
which variations are not heritable and in our theory may be account- 
ed for through the characteristics of b;. We are speaking of heritable 
changes where a discontinuity manifests itself because there are no 
intermediate forms between the unchanged and the few changed. The 
biologist calls these discontinuous changes mutations. We will inter- 
pret these mutations as changes in g which may occur spontane- 
ously in any body cell and changes which the offspring of that body 
cell will possess. These changes are truly spontaneous and intrinsic for, 
as E. Schrédinger (1945) shows, thermal agitation, which is always 
present at the temperature at which living substance is maintained, is 
able to account for these spontaneous mutations. 


The cancer viruses (Rous, 1935-1936) and the ‘milk factor” 
(Bittner, 1942), which is probably of a nature similar to that of can- 
cer viruses, add a degree of complexity. We shall count cancer virus 
as a contributor to g. J. Needham (1942) makes the observation that 
the objection to the virus theory of cancer lies in the fact that it pri- 
marily emphasizes the exogenous origin of uncontrolled specific pro- 
liferation. He wonders whether the difference between transmissable 
and non-transmissable tumors may not be any more mysterious than 
the capacity which some enzyme systems possess and others do not, 
of coming out into cell-free extracts of tissues. Malignancy (Berill, 
1948) may be one of the many properties of a cell induced in host 
cells by its non-cellular extracts. Similar extracts of non-malignant 
cells might have similar induction properties, but without simultane- 
ous conveyance of malignancy and consequent growth they remain 
undetected. There are many properties which the gene and virus 
possess in common and the possibility must exist of the virus being 
a sort of gene-like body which can survive transmission through a 
cell-free extract and which once becoming part of the cellular mech- 
anism remains that (Fischer, 1946). Of interest in this direction is 
that transmission by cell-free filtrates of fowl tumors induced by tar 
has been reported (McIntosh and Selbie, 1939). 

One of the ideas pertinent to neoplasm formation is that a tumor 
represents an escape of cells from a combination of injurious agents 
which would otherwise eventually cause death of the cell. The escape 
is accomplished by the acquisition of a new form of metabolism by 
the cell which cannot be adversely affected by the injurious agents 
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(Moore, 1944). This may be a means by which such a large number 
of cells are born and killed that the probability of a single cell com- 
ing into existence with a mutation bringing a sudden increase in.g 
is vastly increased. This mutant will have abnormal hereditary ten- 
dencies and altered structures. The new metabolic process in the cells 
is neither normal nor subject to the usual restraints, and the growth 
and multiplication of the cells proceeds autonomously. The assump- 
tion of a higher glycolytic coefficient may be one of the alterations and 
one which has been shown (Rashevsky, 1940) adequate to produce 
many of the characteristics of neoplastic growth and cell division. 

Considering g further, we will mention the growth-promoting 
abilities of embryo extract; as well as of tumor cell extracts. Allied 
with it is the observation that rapidly growing embryonic tissues fre- 
quently show pyknotic nuclei and degenerating cells (Needham, 1942). 
It seems that cells have extractable products, which also appear on 
cellular degeneration, and have the power to increase what we call 
the value of g assigned to the cell. 

Our J; is related to the concept of the individuation field in ex- 
perimental embryology. Cancerous growth is an escape from the con- 
trolling pattern or field. 

S. Bayne-Jones et al. (1937) have asked the very important ques- 
tion, “Have they [malignant cells] become ‘fast’ to the conditions 
which normally control cell growth in the body, or is there a break in 
the internal control mechanism of the cell, or is there a loss in body 
control of cell activity? This, the core of the problem, has been al- 
most entirely neglected.”’ We propose to show that all these processes 
together are operative in the production of the cancer cell. 

According to equation (3) positive growth, of a neoplastic type 
or otherwise, will occur only when 


To satisfy equation (41) we must either increase niG; (which is equal 
to n,aig, by definition of G;), or decrease J; , or do both. We will de- 
scribe cases of each of these theoretically predicted methods of car- 
cinogenesis as demonstrated in the experimental laboratory. 

An increase in the n,;a:g9 term is most significantly due to a 
change in g. Such an increase in g has been demonstrated in the neo- 
plastic transformation of normal cells grown in tissue culture. In the 
work of A. Nettleship and W. R. Earle (1943) mouse fibroblasts were 
treated in tissue culture with methylcholanthrene. Cultures of such 
cells, when transplanted subcutaneously into the same strain of mice 
from which they were derived, grew there with all the characteristics 
generally ascribed to cancerous tissue and eventually caused the death 


80 ORGANISM GROWTH AND CANCER PROBLEM 


of their hosts. That the change was in niaig is evidenced by the fact 
that in the tissue into which the cells were transplanted there was 
obviously approximately the same value for J; as in the tissue from 
which the cells were originally extirpated. 

In illustrating the affects of decreased J; we turn to the results 
of A. Fischer (1937), who by transplanting pieces of various tissues 
into abnormal sites in the same animal from which the pieces of tis- 
sue were removed, was able to produce neoplasms from this trans- 
planted tissue. Apparently the inhibition (measured by J;) to the 
proliferation of the cells which was exerted in their normal sites was 
not so great in some of the new locations into which they were trans- 
planted and, therefore, the size of the value of nia:g which remained 
constant was sufficiently great to cause a rapid growth in the new 
sites. 

As a case of an elevation of n;@;g and a depression of J;( = bitin 
[n; — »;] — Gi) occurring simultaneously, reference may be made 
to the work of K. Yamagiwa and K. Ichikawa (1918), in which epi- 
dermoid carcinoma developed on the ears of rabbits painted with coal 
tar. Regarding the effect of tar on J;, it is known that under certain 
conditions methylcholanthrene can be a destroyer of tissue, and can 
cause cell death (Earle and Voegtlin, 1938). With this cell destruc- 
tion there would be a decrease in n; and perhaps also a decrease of 0; , 
with resultant decrease of J;. With random cell death those cells sur- 
viving might absorb some of the diffusable products of their dead 
neighbors, with other changes in g of the cells of the tissue occurring 
by mechanisms described previously as means of elevating g in a cell. 
There is a vast variety of substances which have been shown to be 
carcinogenic and which are certainly very diversified in nature. J. P. 
Greenstein (1947) mentions, too, bakelite discs, zinc salts, x-radia- 
tion, ultra-violet radiation, phthalate buffer in HCl, estrogens, cys- 
ticercus, etc., as carcinogenic agents. The simplest generalization 
which might be drawn from this is that any agent which may be ap- 
plied over a prolonged period of time and which causes a selective 
necrotizing process of cells but spares some cells may cause those 
cells surviving to undergo malignant change. In this connection it 
might be interesting to mention the action of benzol on bone marrow 
(Goodman and Gilman, 1941). Benzol, because of its destructive ef- 
fects on bone marrow cells, has previously been used in the treatment 
of the lukemias (tumors of bone marrow). Paradoxically, chronic 
benzol poisoning has produced lukemia as an end result. Likewise, 
x-rays are one of the accepted treatments for cancer, yet it is well 
known that exposure to x-radiation can produce cancer in a previously 
healthy tissue. This same idea may be extended to many other car- 
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cinogenic agents. 

The influence of age on the rate of growth of a tumor is well 
known, with tumors in early life tending to progress much more rap- 
idly than do the same tumors in later life. Assuming all other quan- 
tities are the same, this may readily be seen from: 


dn;/dt = niG; — b,tn(n; — vi) my viGi ) 


which gives the rate of growth of a tissue region (dn;/dt), no matter 
what be the cause of the growth, as a function of age (t). This equation 
under other conditions can predict the falling off of the wound-healing 
rate with increasing age. 

The incidence of cancer seems to be an increasing function of age 
despite the inaccuracies inherent in our present statistics. Following 
the approach of N. Rashevsky (1945), we will show how the theory 
may account for this increased probability of occurrence. Using equa- 
tions (41) and (5) we see that for growth to occur at all we must have 


niGi = [b;in(n; — vi) Se viGi] f 
It follows that 
niG; + 1G; > bitn(n, — »:), 


and 


G; N; == Vj 
Be < ( ) (42) 


in dd reams Ue 


It is known by the statistical treatment of nearly all fluctuating 
characters that among the members of any large group of organisms 
we obtain the form of the normal curve of probability (Woodruff, 
1941). Of the different quantities used in our equations that vary 
significantly among the different members of a species, we have g and 
b;. At present we might suppose that g varies in the discontinuous 
way which hereditable variations do. Quantity b; will be considered 
as one of those fluctuating continuous variations which has an aver- 
age value in the species of 6; and where the probability of b; in a par- 
ticular individual deviating from 6; by an amount (6; — bi) is given 
by 

N (0; — bi) = Ae-a(bi-bi)? , (43) 


We get the probability Wi(t) of a cancer occurring in the 7th 
tissue region of an individual at a given age from the following equa- 
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tion by substituting equation (42) into (43), giving us 


Gj NiEtVi 
W;(t) —yerk- 3] : (44) 

It can be seen that W;(¢) will increase with age, since 
Gi (= = Vi ) 
tN \ 2; — 9; | 
is an increasing function of age. In the numerator we only have n; 
decreasing with age, while in the denominator we have (7; 2 yi) get- 
ting very small with age, with a decreasing n and a ¢ increasing rela- 
tively slowly compared to the decrease of (ni — ¥:). 

Grateful acknowledgment is made to Professor N. Rashevsky for 
his inspiration and invaluable counsel, and also to Dr. Anatol Rapo- 


port for his helpful suggestions and criticisms during the preparation 
of this paper. 
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ADDENDUM 


According to W. R. Earle (personal communication dated April 
2, 1948), it has been possible to grow the single cell of one cell type in 
tissue culture, with the difficulty in accomplishing this probably due 
to the unsatisfactory nature of the culture medium used in tissue 
culture. 

This recent accomplishment would not seem to make untenable 
our concept of v; as the minimum number of cells of type 4 which 
must be together in aggregation when any cell division will occur at 
all, insofar as the intact organism is concerned. Growth rate in tissue 
culture is usually more rapid than that seen in the organism in vivo. 
The body fluids would then not seem as satisfactory (for furnishing 
media for positive cell growth) as some of the artificial fluids used. 
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In previous tissue culture growth experiments, single cells could not 
be grown, with a minimum number of cells required for growth to 
occur at all. Reasoning from these facts, it would seem that in vivo a 
larger group of tissue cells would be required to initiate a positive ; 
growth change. See 
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CAN ELASTIC STRESSES IN GELS CAUSE THE ELONGATION 
AND DIVISION OF A CELL? 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


Some metabolites present in the cell may affect its colloidal structure. 
If the cytoplasm is in the state of a gel, elastic stresses may appear as 
a result of non-uniform concentration of some metabolites which affect 
mechanical properties of gels. This paper investigates the question 
whether mechanical stresses, thus produced, can result in such a defor- 
mation of the cell which may lead to a division. The result of this study 
indicates that this is not possible. Elastic stresses, produced by non-uni- 
formities of concentration of metabolites which affect the colloidal struc- 
tures of the cytoplasm, may produce only a general dilation or contrac- 
tion of the system, without change in shape. 


The success of the diffusion forces theory of cell division (Ra- 
shevsky, 1940, 1948; Landahl, 1942) in predicting with accuracy the 
time course of the elongation and constriction of freely dividing de- 
membranated Arbacia eggs (Buchsbaum and Williamson, 1943) has 
led to a certain neglect in the study of other possible causes of cell di- 
vision. The only other quantitatively treated possibilities were those 
of variable surface tension (Rashevsky, 1931, 1938; Shimbel, 1946), 
and the purely formal approach based on the principle of maximum 
energy exchange (Rashevsky, 1943; Stanton, 1944). 

While these approaches also lead to a mathematical description 
of the dynamics of elongation, they do not directly give us any indi- 
cation as to the order of magnitude of the quantities involved. In 
that respect the diffusion force theory has a definite advantage. (Ra- 
shevsky, 1938; 1940). 

Because the diffusion force theory gives us only the description 
of the overall process of the division of a cell as a whole, and is in- 
adequate for the description of the phenomena of nuclear division 
(Rashevsky, 1938; Schrader, 1944), all other possibilities should be 
studied mathematically in a systematic way. 

The protoplasm is a complex colloidal system which may under- 
go various sol-gel transformations. Substances dissolved in'a gel] in 
general may cause phenomena of hydration, or other phenomena 
which result in the change of the specific volume of the gel. Non- 
uniformities of the concentrations of a dissolved substance may thus 
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produce non-uniformities in expansion or contraction of a gel, thus 
causing mechanical stresses. The question arises whether such me- 
chanical stresses may produce an overall elongation of a cell with a 
resulting division. The question becomes of special interest because 
the biological importance of mechanical stresses in gels has been em- 
phasized in the field of embryology (Weiss, 1934; 1940). 

The presence of protoplasmic streamings during cell division 
would seem to indicate its fluidity and any assumption of stresses 
may appear absurd. The question is nevertheless worth studying, 
because the protoplasm is known to sometimes possess both the prop- 
erties of an elastic solid, and of a viscous liquid. 

In developing a mathematical theory of such a contradictory 
physicochemical system, we will have to use somewhat contradictory 
methods. In computing the stresses, we shall consider the system 
as a solid. From the stresses, we shall compute the volume and sur- 
face forces acting on the system. Then we shall assume that under 
the prolonged effect of such forces, the system behaves as a plastic 
body, and thus compute its rate of deformation. 

Let the specific volume of an element of our system be V, in the 
absence of any dissolved substances. Let a substance be present 
in non-uniform concentration c(x, y, z), and let that substance 
cause a dilation of the system. If a denotes a constant and ¢(c) any 
arbitrary continuous function of c, then the linear expansion may 
be represented by 1 + a ¢(c), and the volume V, changes to 
Vo[1 + 3 a ¢(c)]. The problem of computing the stresses due to 
such an expression is formally identical with the problem of comput- 
ing thermal stresses in a solid with a coefficient of linear thermal 
expansion « in a non-uniform field of temperature ¢. Denote by E 
the modulus of elasticity and by o the Poisson ratio. The elastic 
stresses and strains in such a system are the same (Timoshenko, 
1934) as produced by a system of volume forces F: 


F =— grad Ane. =—orad Oso Geass ane : (1) 
per unit volume, and by a system of normal surface forces 
ae eae (2) 
1—2c 


the outward direction of the normal being taken as positive. 
We now ask the question: if the forces (1) and (2) affect a plas- 


tic body of viscosity » what will the average rate of plastic flow in 
a given direction be? 
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To answer that question, we may use Betti’s theorem, as has been 
done before (Rashevsky, 1940; 1948). It must, however, be noted, 
that whereas before we considered the plastic body as incompressible, 
and therefore set o = 4, we cannot do that now. In fact, putting 
o = ¥ will make both forces (1) and (2) infinite, which is rather 
obvious physically. Hence, we must use Betti’s theorem in its most 
general form. Denoting the overall length of the system in the z-di- 
rection by 1, , by V—the total volume of the system—by X , Y, Z, and 
X,, Y,, Zy the components of F and F, respectively, we now have: 


1dl, 
ie Sree 


(3) 
2 ae + 2X,)] dS 


\ 


Even though F = grad @, the volume integral now cannot be trans- 


formed completely into a surface integral, as is the case for o = 4 


(Young, 1939). 
Integrating by part the first term of the volume integral in (3) 
we find now 


ff f-=a- snieeet Ge) RES 


Similar expressions are obtained for the other two terms. The whole 
volume integral now becomes equal to 


| foc cos(y,z) —o[ycos(y, y) 
s (5) 


+2zcos(y,2)]}.dS— (1—2<¢) JJ Jee. 


From equation (2) it follows that 
X,=F, cos(v, x) =— ® cos(r, x), ete. (6) 


Introducing expressions (5) and (6) into equation (3) we find 
that the two surface integrals cancel each other out and thus, because 
of equation a we obtain: 


et oy av JJ cea) ete 
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Exactly the same expressions are obtained for (1/l,) (dl./dt) and 
(1/l,) (dl,/dt). Hence, the shape of the body remains unchanged. For 
¢ > 0. an average dilation occurs, for ¢ < 0—an average compression. 
No overall change in shape, such as an elongation and constriction 
will occur in this case. 

It may be argued that this negative result is due to the contra- 
dictory assumptions made about the simultaneous behavior of the 
system as an elastic solid and a plastic fluid. But the negative result 
holds a fortiori if we make more rigorous assumptions. If the sys- 
tem behaves only as an elastic body, no plastic flow is possible. By 
now applying Betti’s theorem to calculate the average elastic strains, 
we find again that only an average constant dilation or compression 
remains. If, on the other hand, the system is a purely plastic body, 
no stresses are possible, and hence no forces due to stresses. 

It must be kept in mind, however, that Betti’s theorem gives only 
the average rate of plastic. deformation in a given direction, or, for 
an elastic system, only the average strain in that direction. From the 
fact that these average values are zero or the same in all directions, 
it does not follow that the same theorem holds for every element of 
volume of the system. The system of forces (1) and (2) may in gen- 
eral produce internal plastic flows, leaving the overall shape of the 
body invariant. The calculation of these local flows is, however, a 
problem of an entirely different order of mathematical difficulty. It 
is quite possible that such internal local displacements may play a 
role in such cases of mitosis which leave the cell as a whole unde- 
formed. 
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CONDITIONING STIMULUS UPON THE RESPONSE 
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The strength of a test stimulus, which is just adequate to produce 
a response after a sub-threshold conditioning stimulus, is calculated on 
the basis of Rashevsky’s two-factor theory as a function of the strength 
and duration of the conditioning stimulus. The results are compared with 
available data for one of a family of curves and found to be in satisfac- 
tory agreement. 


If a nerve is subjected to a sub-threshold stimulation, its response 
to a test stimulation is altered. In the present paper we shall con- 
sider, in terms of the two-factor theory of N. Rashevsky (1938), the 
effect of a conditioning stimulus on a cathodal test stimulation when 
both are constant currents. This problem has been investigated ex- 
perimentally by J. Erlanger and H. 8. Gasser (1937). 

At the time t = 0 let a constant sub-threshold current f.R , R be- 
ing the rheobase, be applied and maintained for a time;. Att™=v7z 
let the current be suddenly changed to (f. + A)R so that A is the 
fraction of the rheobasic current added at t = 7 and maintained in- 
definitely. We shall consider such values of A as are just adequate to 
initiate excitation at some time after t — 7. Let excitation occur at 
atimet=t# +7. 

Following N. Rashevsky (1938), we may write the differential 
equations governing the excitatory state « and the inhibitory state 
j in terms of the constants K, M, k and m and the initial values e 


and Jo: 


T= KIR—kle— a) ’ (1) 
4 Mp R— mj — 50) (2) 
dt c. OJ 


Since the current in the interval between zero and ; is constant, 
we may write for «(7) 
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Kf 
e(r) =& + u 


(1 — e*7) (3) 


and a similar expression for 7(r). For any time t > + the current is 
constant, and equal to (f, + A)R. Hence, since « = e(r) is given 
by equation (3), we may write, if ¢ =t—7, 


ole 
e(t > 7) =& ual (1 —e*")e*" 
| K(f,+ A)R a 
~ pe eee (1—e*"), 


A similar expression holds for 7(¢ > 7). 

Excitation occurs at the time t = ¢* when e = 7. Hence setting 
e(t > +) of equation (4) equal to the similar expression for 7(¢ > 7) 
and then introducing t’ = ¢*, we obtain an expression relating 7, t’, 
f. and A: 


. _ KTR 
6G tO) IG tO) =a mie tee ha) ee 
Mf.R c 
_ Mick er Sees. Emr) em + K ee F (1 aes re (5) 
M(f.-+A 
_ Mich R Sara 


Since A is by definition the smallest value necessary for the test stimu- 
lus to produce excitation, t* must also determine a maximum for « — 7. 
Thus for t’ = ¢*, d[e(r + t) —7(r + t')]/dt =0. Hence we may 
solve for t* and obtain 


1 K\ fA+f 6% 
log ( = ) =.) ; (6) 
k—m M/\A+fe™ 
The cathodal rheobase, R , may be determined from the equation 
(Rashevsky, 1938) : 


Jone fy 7K. M \¥/-m) M M \/ (km) 
EG -EE UY] 


Substituting values for ¢* and R from equations (6) and (7) into 
equation (5) we have 


Cee 
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Mk 
(mn — = ) re! —fe— A) (A aie fice *™) / em) 


Mk: M m/(k-m) M kK/(k-m) ; 
ita ee) (8) 
xX {(A + J e@-™7 )*/ em) = (A + Teen | ; 


which determines the relation between A and the independent vari- 
ables f, and z in terms of three constants (M/K), k and m. If we 
assume that K/k = M/m (normal accommodation), equation (8) re- 
duces to ; 

(A+ f.e"")* = (A + fe)”, (9) 


so that A(f., 7) is determined by the two constants k and m. 
In Figure 1, 4 is plotted against ; for f. = —0.5, 0.5 and 0.9, 


--- k=1180 

m= 59 

© pata FROM 
ERLANGER & GASSER 
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FIGURE 1 


using equation (9). The dotted lines are obtained by. using k = 1160 
and k/m = 20. The solid line is given to illustrate the dependence 
upon k/m. In this case k = 750 and k/m = 10. The points are ex- 
perimental values by J. Erlanger and H. S. Gasser (1937), Figure 44, 
page 95. The data are represented particularly well by the solid curve. 
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If k and m are determined independently from intensity-time data 
with constant current and from data with alternating current, then 
the curves are completely determined unless K/k # M/m. 

When K/k # M/m, equation (8) must be used. From this equa- 
tion it-is not easy to determine A for given f, and 7 even if the con- 
stants are known. Define p and o. by the equations: p = k/m and 
M/m = (1+ a)K/k. Ifp< 0.1, a< 0.1, and 4 ss 0.1 (which holds 
at least if f, < 0.9) , then, since a? ‘is approximately equal to 1 + p 
log a, we may y write very roughly: 

A=1—f.(e"" —a)/(1—a) ++ p(1 +a) log(4 + f.e*).. (10) 


We may take as a first approximation A, to A the expression giv- 
en by the right-hand side of equation (10) but with the third term 
omitted. Then substituting this value into the right-hand side of equa- 
tion (10) we have the second approximation A, to 4 

Coo 


A,=1—f;, 


—— ty, 
(11) 
a 
Bey Pose bt 


em — 
+ p(1 +a) log (1—f-. ji 


This expression is actually very satisfactory if p < 0.1, a < 0.1, 
f- < 0.9 and if A < 0.8. The last restriction may be removed if A, is 


~ 


substituted for A on the right-hand side of (10) and a third approxi- 
mation is obtained. Further substitutions do not help since equation 
(10) is an approximation to (8). The values obtained in this man- 
ner are almost indistinguishable from those given in Figure 1. 

The parameters may be estimated from experimental data as fol- 
lows: When 7; becomes very large we have from equation (8) 


A 5s) et oa elo 
| (3a) eS Saag 


afc es 


+, 
1+ plogp—p 
Thus a may be estimated to a first approximation by neglecting the 


terms in p. When - is large enough (m+ > De eqususn (11) reduces 
to 


A,= A, —f.(1 =p) er (Lay, (13) 


where A., is A2(+ = «) and depends upon f, . Introducing the experi- 
mental values for A. and f, and the estimate of a, one can estimate m 
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to a first approximation by setting p = 0 in equation (18). With these 
values of a and m, k can be roughly estimated by trial and error using 
equation (11). With this estimate, and since p = m/k, both a and m 
can be adjusted. Then with the new values for a and m, a corrected 
value of k can be obtained and A(f,, 7) can be computed from (11) 
or (8). 
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A MECHANISM FOR OPTIC NERVE CONDUCTION AND 
FORM PERCEPTION: II 


JAMES T. CULBERTSON 
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This paper continues the analysis of a hypothetical device for “bottle- 
neck” optic nerve conduction and form perception. The neural mechanism 
described in this paper provides for a form perception which is invari- 
ant for transformations of size, translation and rotation on the retina. 
The number of neurons at different levels in the mechanism is considered. 


I. Form perception invariant for changes in size, position, ori- 
entation of the retinal image. 


In a previous paper (Culbertson, 1948, hereinafter denoted by 
MI) we began the description of a hypothetical retina and nerve net. 
This description is continued in the present paper, which pre-supposes 
a knowledge of MI. In MI we showed that if any pattern (activity 
matrix) occurs in 2’ (any one of the »y = 100 square arrays of 225 
neurons each at level Z) at any time ¢, then the same pattern oc- 
curred at t — 1 in some set s* of neurons in the synchronized retinas 
D or D’. The set s* consists of all the neurons within some one of the 
6.4 X 107 mappings on D and D’, the mapping subdividing s* by small 
squares into 225 subsets whose firings at t — 1 collectively determine 
the activity matrix P(s* , £— 1), which is called the pattern in s* at 
t — 1*,i.e., for any values of tand]—1,2,---» we have P(z', t) = 
P(s*,t —1) always true for some value of «. 

Remembering that all the y sets at level Z are identical in struc- 
ture, and the » sets at level H are identical in structure and also in 
the contacts they make with the Z level, we have the following re- 
sult. Let any light image be projected on the retina so as to cause 

* A few of the non-upright mappings of the 4 smallest sizes, (those with 
small squares measuring 1, 1.1, 1.21, and 1.331 units on each side) and a few of 
the large non-upright mappings at the periphery, will have some small squares 


not containing any neurons, so that in these cases s* contains some null subsets. 
In the pattern P(sk, t — 1), the activity value of any null subset s;; is defined 
as zero. Considering the total number of mappings evaluated in MI, we see that 
a negligible number of small mappings gives identical neuron sets. Thus, rotat- 
ing any smallest upright mapping through 3.6° or multiplying its area by 1.1, we 
find the same neurons enclosed, and similarly for some other special cases involv- 
ing very small mappings. If the same neuron set is given by two mappings, the 
mechanism funnels the set twice during each cycle period, once for eaeh mapping. 
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some pattern in z’ at ¢ such that P(z', t) comes from an s“ having 
no null subsets (see previous footnote). Also, let all neurons in oe 
which fire at ¢, say there are w of these, have an excitatory end bulb 
on some neuron h, at level H and all not firing have an absolutely in- 
hibitory bulb in h.. Thus, by the structural identity mentioned, there 
are y neurons similar to h., one above each 2’. Each is designated 
h,. Let each h, have one bulb on some neuron f., whose threshold is 
1. Then there are N- = 6.4 X 107 changes in position, size, orienta- 
tion of the image which likewise cause f, to fire if each projection 
lasts T = mké = 1280 sec.* Here there are no other end bulbs on h, 
or f., the threshold of h, equals the number of excitatory bulbs on it, 
and NX, slightly less than N, is the number of distinct mapped sets 
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s* without null subsets. As our example, let us again take pattern 
P, and define P; as the activity matrix in Figure 1. Thus 


fi(é+ 2) =(oh(t + 1)) = [P(2',t) = Pr, 


where 1=1,2,---». As our “T-shaped” image we take a solid T 11c 
units across, 11c units high, and 3c units thick through the horizontal 
and vertical bars so that its area is we? where \/2 =c = 3334. In Figure 
1 this is shown superimposed on an upright mapping (see the next 
to the last paragraph in MI). The unit of distance, as mentioned pre- 
viously, is the distance between two consecutive neurons in any row 

or column of A or D. % 
For any T-image projection (thus defined) causing f, to re, we 
notice that small variations of its position, size, orientation likewise 
* The necessary projection time and hence also the reaction time is, of course, 


much too large as mentioned in MI. This will be reduced by an accessory “atten- 
tion” apparatus involving a small number of neurons. 
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cause f. to fire and that, because of the small gradations in the map- 
pings, it is possible to move the image continuously from one position, 
size, orientation causing f, to fire to any other causing f, to fire (if 
c¢ > \V/2) so that every intermediate variation of position, size, orien- 
tation can cause f, to fire. 

As explained in MI, there were mappings of 123 different sizes. 
Their areas were 225(1.1)¥ where y = 0,1, 2,--- 122. The center 
squares of equal sized upright mappings did not overlap but filled a 
square area whose center was the center of D. The number, Ny, of 
upright mappings on D of area 225(1.1)’ was a maximum such that 


R=N,= (10+)? 


where FR was the greatest number of upright mappings of area 
225 (1.1)” that could be mapped on D without any overlapping of their 
center squares. For y S 81 the value of Ny was restricted by the right- 
hand member of the above expression; for y > 81 it was restricted by 
the left-hand member. Hence, all upright mappings on D of area 
225 (1.1)7 are contained within an upright restricting square Uy which 
is (24 + y)?(1.1)’ units in area centered on the center of the retina, 
so that mappings less than 225(1.1)® in area are restricted by these 
squares to less than the entire retina, those greater are restricted by 
the boundaries of the retina. Thus, the minimum size of a T-image 
which can cause h, to fire is a function of retinal position. The small- 
est T-image which can be translated from any projection in U, effec- 
tive in firing f, to any other effective projection in U, continuously, 
so that each intermediate projection can be effective, has the area w. 
In general, a T-image of area w(1.1)7 can be translated from any ef- 
fective projection in Uy, to any other effective projection in Uy con- 
tinuously, so that each intermediate projection is effective. Also any 
T-image of area A = w(1.1)’, where y > 8, projected within the reti- 
nal surface and within Uy for T = mké , will cause f. to fire. That is, 
any T-image (as defined above): projected on the retina for mké 
in any position, size, orientation causes f. to fire if its area 
A = ow(1.1)” and it falls entirely within the retina and within U,. 
This shows a theoretically possible mechanism which indicates in part 
a solution to the gestalt problem of invariant response to an image 
of a given shape under position, size and orientation variations. The 
T-image is, of course, just. an example. We must indicate how varia- 
tion (within any limits we wish to assign) in the shape of the image 
may result in invariant response. This is done by relaxing the firing 
conditions for h,. Also, we must indicate how there are many differ- 
ent possible classifications of patterns into forms. A form is defined 
as any set of patterns firing the same neuron at level F' (the term 
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“shape” will refer to a property of the image). The number u of 
forms which can occur in the mechanism equals the number of h neu- 
rons contacting each z’ which, because of the way levels F and H are 
connected, as already described, also equals the number of neurons in 
level F’. 

It seems to be, to some extent, a question of how far we want the 
shape as well as size, position, orientation of each image to be alter- 
able and yet produce an invariant response, and what forms, and how 
many, we want this mechanism to respond to invariantly. If we de- 
cide what we want in regard to these questions, then, within obvious 
limitations, we can get the desired results. There are limitations on 
the number of neurons, the number of dendrites and axons per neu- 
ron, and the number of synapses per dendrite or cell body. 

Arbitrarily we assign u = 50,000,000 forms to the mechanism, 
that is, 50 million classes of patterns in z’ such that the same neuron 
f is fired by all patterns within any one class and is not fired by any 
patterns not in that class. (Certain of these classes have some pat- 
terns in common, as will be indicated.) The order of P(z', t) fixes 
the number of pattérns in this mechanism at 2°”. 


II. Neuron economy versus dendrite, axon, and geometric econ- 
omy. 


For simplicity in explaining the principles involved, we placed 
no limit on the number of dendrites and axons a neuron could have, 
nor on the maximum distances between dendrite (or axon) endings 
of the same nerve cell. Economizing in this nerve net to a maximum 
of approximately 200 dendrites, to the same maximum for axons, and 
to plausible limits on the distances mentioned, the consequent altera- 
tions involve an increase in the number of neurons in the net (but 
otherwise no important changes). It is necessary to show that this 
economy still allows a neuron economy giving a reasonable number of 
nerve cells. 

In Figure 1 in MI, the principle layers in the nerve net are shown. 
The elements up to level C are in the retina and optic tract. Let us 
consider them first. 

There are w = 2.5 X 10° receptor neurons (level A) and these 
are divided into m, = 25 subsets, each subset A* containing n, = 10° 
neurons. The cycle contains m,k = 100 neurons, so that if there is 
just this one cycle (call it J,), as assumed for simplicity in MI, then 
each of the m, neurons designated x, must have 10* end bulbs con- 
tacting level B , which is not plausible. Instead of this, we will con- 
struct the net so that there are 10,000 similar cycles, J, , Jo, -- J 566 
reducing the number of these end bulbs for x to 100, which ives 
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cycle neurons. The cycles are synchronized by allowing any one mem- 
ber of J; to have an end bulb on its successor in Jj., or they may be 
connected in parallel. 

Although the receptor neurons may in theory be assigned in any 
way to the m2, subsets, some ways of assigning them give rise to ana- 
tomical improbabilities. For instance, suppose that we divided the 
retina into m, equal square areas, each containing ~) neurons and 
assigned all neurons within any one square to the same subset A’. 
This requires each cycle to cover widely separated points throughout 
almost the whole retina as it must extend to each square, and this 
seems improbable. Also, each cell at the B level must have its m, den- 
drites extend to widely separated points as each must enter a differ- 
ent square, and this likewise seems improbable. However, we may 
divide the retina into m square areas each containing m, neurons and 
assign each neuron within any square to a different subset.. Now each 
cell at the B level has its m, dendrites extend only to points within a 
square area 1/10° the area of the retina and each cycle may be con- 
fined to an area only 1/10‘ the area of the retina. It is clear that for 
any hypothetical nerve net such geometrical economy must be satis- 
fied, as well as the numerical economies. 

There are % + ™,. = 1,000,025 = 10° neurons in the optic tract 
of the mechanism (see level VX) so that in the retina and optic tract 
there are approximately 2.7 < 10% neurons. 


rete 
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Element c;* has the same relative position in C as @;* has in A, 
so there is geometric as well as numerical economy in the contacts 
made by the end bulbs of v;. In MI, however, for simplicity in ex- 
plaining the principles involved, we gave each x* contact with 2 x 10° 
cells, one excitatory on each member of C and one inhibitory on each 
circuit between C* and D*, which is neither geometrically nor nu- 
merically economical. We could achieve these economies by introduc- 
ing a set of internuncials after each x‘, but there is a simpler sym- 
metrical method shown in Figure 2. 

We may remove the entire neuron set X (x1, x?.--- g--- x) 
from the nerve net, together with that part of x, leading only to 2* 
(see Figure 1 in MI) and insert a set of cycles Jy’, Js’ --- J’10,000 €X- 
actly similar to the set J,, Jo --+ J10,000 and associated with level C . 
The neuron in J,’ corresponding to £ in J, is designated f’ and 6 ex- 
tends through the optic nerve to have an end bulb on f’, but there is 
no other connection between the two sets of cycles. 

At any time t (> 0) there is one and only one cell firing in each 
of the 2 X 104 cycles. The neuron in each J’ cycle corresponding to 
x, in each J cycle is designated x*. All x, fire at the same time and 
all x* fire at the same time (1), but the latter fire one 6 after the 
former (2). 


IT 2) (ft). = 2. (6)> also. IL a6 = (1) 
v(t) =x (t—1). (2) 


A trivial adjustment necessary in the above is to take any time that 
all 8 fire as the time origin, i.e., replace Proposition 1 in MI by ZI £ 
(0). Then all other numbered propositions in MI remain true. The 
set J’ has the same relation to level C as J has to B, so that if x in 
J, has an end bulb in b,*, then x in J,’ has an end bulb on ¢;*, but in - 
addition x* has an end bulb inhibiting the circuit above ¢;>. 

This adds 10° neurons associated with set C and the latter con- 
tains w = 2.5 X 107 neurons while there is a negligible decrease of 
24 fibers in the optic nerve. 

In a subsequent paper we (1) complete the economy analysis 
from level C through level F', (2) indicate the classification of the 
2° patterns into u = 5 X 107 forms, and (3) construct the “atten- 
tion” apparatus to reduce the mechanism’s reaction time. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 
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